UMD-PP-10-020 
arXiv:1011.5090 



Current and Future Constraints on Dark Matter from Prompt and Inverse-Compton 
Photon Emission in the Isotropic Diffuse Gamma-ray Background 

Kevork N. Abazajian 1 ' 2 ^] Steve Blanchet 2 ' 3 [j] and J. Patrick Harding 1 ' 2 ^] 

1 Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697 USA 
2 Maryland Center for Fundamental Physics & Joint Space-Science Institute, 
Department of Physics, University of Maryland, College Park, Maryland 20742 USA and 
3 Instituto de Fisica Teorica, IFT-U AM/ CSIC Nicolas Cabrera 15, UAM Cantoblanco, 28049 Madrid, Spam 

(Dated: February 9, 2012) 

We perform a detailed examination of current constraints on annihilating and decaying dark 
matter models from both prompt and inverse-Compton emission photons, including both model- 
dependent and model-independent bounds. We also show that the observed isotropic diffuse gamma- 
ray background (DGRB), which provides one of the most conservative constraints on models of 
annihilating weak-scale dark matter particles, may enhance its sensitivity by a factor of ~ 2 to 3 
(95% C.L.) as the Fermi-LAT experiment resolves DGRB contributing blazar sources with five years 
of observation. For our forecasts, we employ the results of constraints to the luminosity-dependent 
density evolution plus blazar spectral energy distribution sequence model, which is constrained by 
the DGRB and blazar source count distribution function. 

PACS numbers: 95.35. +d,95.55.Ka 



I. INTRODUCTION 

The existence of cosmological dark matter is well estab- 
lished by observations of galaxy clusters, galaxy rotation 
curves, the cosmic microwave background (CMB) and 
large-scale structure, though its nature remains a funda- 
mental problem in cosmology and particle physics. There 
exists an abundance of particle candidates which could 
account for the dark matter (for a review, see, e.g. [1]). 
For a class of particles with weak-scale interaction and 
weak-scale particle mass, weakly interacting massive par- 
ticles ( WIMPs) , their production in the early Universe in 
thermal processes naturally produces the observed dark 
matter density, largely independent of the particle mass. 

Thermal freeze-out predicts a canonical annihilation 
rate of (<tai>) ~ 3x 10 -26 cm 3 s _1 . This predicted annihi- 
lation rate in standard model channels leads to energetic 
gamma-ray production in the hadronization of quarks, 
the Higgs or gauge bosons, through bremsstrahlung in 
the case of the lighter leptons, or directly to two gammas 
through higher order processes. The diffuse gamma-ray 
background (DGRB) was forecast to be one of the most 
robust constraints on annihilating WIMP dark matter 
[2]. Because of a more conservative model for the ex- 
tragalactic dark matter signal, the conservative limits on 
dark matter annihilation presented by the Fermi-LAT 
Collaboration from the 1-year observation of the DGRB 
[31 0] were weaker than prelaunch estimates [2] ■ On the 
other hand, it was shown that the DGRB has an irre- 
ducible contribution from the Milky Way Galactic dark 



* kevork@uci.edu 

t steve. blanchet@epfl.ch 

* [EaJcTD 923@>um d.edu| 



matter halo [5] that is greater in amplitude than the con- 
servative estimates of the extragalactic contribution, and 
correspondingly has more stringent limits over many an- 
nihilation channels [BJ. This irreducible, isotropic compo- 
nent is due to the fact that annihilating or decaying dark 
matter in the Milky Way halo has an isotropic compo- 
nent equal to the minimum of the annihilation or decay 
signal. This minimum is equal to the amount toward the 
anti-Galactic Center. 

The isotropic DGRB has several potential astrophys- 
ical source contributions, including blazars [7HT0], star- 
burst galaxies [TT] and millisecond pulsars j.2] . The only 
model that successfully predicts the shape and ampli- 
tude of the DGRB over all energies is the luminosity- 
dependent density evolution (LDDE) blazar spectral 
energy density (SED) sequence model with an active 
galactic nuclei (AGN) contribution [MlQj]. The SED- 
sequence model matches the shape of the observed blazar 
SED's luminosity dependence [TS]. Reference [T3] re- 
produces well the DGRB as observed by Fermi-LAT, 
while several other analyses under-produce the DGRB 
from blazars. Reference [T3] estimates that > 98% of 
the blazar flux contributing to the DGRB will be re- 
solved in the 5-year Fermi-LAT survey. Prior work that 
under-produces the DGRB uses a single power-law for 
the spectrum of all blazars instead of the observed SED 
sequence for blazars, e.g [16] . Other recent work with 
varied blazar population models, including spectral shape 
variation |17j . possible point source confusion |18j . and 
BL Lac dominance of the unresolved portion [TH] also 
find that a substantial portion of the DGRB could arise 
from the blazar population. 

There may also be an unmodeled, unreduced isotropic 
Galactic component to the DGRB [3j H]. It should be 
noted that two things could happen if there is a presently 
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unremoved Galactic isotropic diffuse component: one, it 
is detected, modeled, and removed, which will make fu- 
ture constraints stronger; or, two, it remains a systematic 
diffuse background, which means our starting assump- 
tion of an LDDE SED-sequence model is not the correct 
model for the DGRB. However, this uncertainty cannot 
be removed without further observational analysis. 

Below, we calculate and show the current constraints 
from the DGRB on dark matter annihilation and decay 
gamma-rays from the prompt as well as inverse- Compton 
components. In addition, when adopting the LDDE 
plus SED-sequence model forecasts of the Fermi-LAT re- 
solved DGRB, future observations will extend the reach 
of Fermi-LAT sensitivity to dark matter typically by a 
factor of 2 to 3. We explore in detail the forecasts on 
standard WIMP dark matter and leptonic-channel mo- 
tivated models [3D] including Asymmetric Dark Matter 
models [21], nonthermal winolike dark matter [22], and 
decaying dark matter [22], as well as models of light (~8 
GeV) WIMP dark matter in scalar dark matter mod- 
els [23J [2Ej- The resolution of the DGRB into point source 
blazars will reduce the DGRB amplitude and ultimately 
enhance the limits on annihilation in WIMP dark matter 
models. 

In a companion paper, Ref. 14 (ABH2), we con- 
strain our adopted SED-sequence model using the ob- 
served DGRB spectrum as well as the observed blazar 
source count distribution function, dN/dF. In agreement 
with Ref. [13], ABH2 found that > 95% of the flux from 
blazars will be resolved with 5 years of Fermi-LAT obser- 
vation. The resolution of the DGRB is in fact similar to 
the resolution of the cosmic x-ray background observed 
by Chandra, which in turn provided stringent constraints 
on decaying light sterile neutrino dark matter |26j . 

The work presented here advances previous work on 
prompt gamma-ray emission in annihilating dark mat- 
ter (e.g., [6]) by including the enhanced constraints and 
sensitivity from inverse-Compton (IC) emission present 
in the Fermi-LAT observation of the DGRB as well as 
forecasts of the improvement of this sensitivity. In addi- 
tion, we go beyond previous analyses of IC emission en- 
hancement of the extragalactic and Galactic signals (e.g., 
[221 [57]) by applying the IC enhancement in the Fermi- 
LAT observation of the DGRB and its forecast improve- 
ment. 



II. THE BLAZAR POPULATION AND 
SED-SEQUENCE MODEL 

The SED-sequence model specifies the cosmological 
blazar spatial distribution and spectrum for a given 
blazar luminosity. It is based on the observed evolu- 
tion of the peak flux in synchrotron and IC emission 
with luminosity. The luminosity-dependent density evo- 
lution model specifics the gamma-ray luminosity function 
of blazars through a fraction of the total AGN popula- 
tion and its x-ray luminosity function. Our blazar pop- 



ulation and SED-sequence model in ABH2 successfully 
reproduces the observed DGRB and blazar source count 
dN/dF. 

Our model is a modification of that by Inoue and 
Totani [5J, and is detailed in ABH2. We provide a sum- 
mary here. The bolometric blazar jet luminosity P and 
disk x-ray luminosity Lx are related by P = lO q L>x- The 
blazar gamma-ray redshift-dependent luminosity func- 
tion is given as a fraction k of the AGN x-ray luminosity 
function (XLF), p 1 {L 1 ,z) = n(dLx /dL 7 )px(Lx , z). We 
adopt the AGN XLF of Ueda et al. [25] • The main fit pa- 
rameter in the XLF is the faint-end slope index, 71 . The 
model also includes a nonblazar AGN component which 
dominates at lower energies, £L < 1 GeV. 

In ABH2, we constrain the blazar population model 
by simultaneously fitting the DGRB spectrum as well as 
the blazar flux source count distribution function dN / dF 
observed by Fermi-LAT [21 [TB]. The best fit parame- 
ters we find are q = 4.191^, log^/lO" 6 ) = 0.38±g;^, 
and 71 = 1.511q o9- These are consistent with previous 
work [5J , though more constrained because we are also fit- 
ting the source count distribution function dN/dF. The 
model reproduces the DGRB and blazar dN/dF, with a 
reduced x 2 / D OF = 0-63. 

Using the dN/dF estimated from a power- law blazar 
spectrum model is not perfect, since the efficiency de- 
pends on this model [T|>]. However, Ref. jT5] tested the 
dN/dF estimate with a non-power-law fit to the blazar 
spectra and found it did not appreciably change the esti- 
mates of dN/dF, adding a systematic uncertainty of 10%. 
We also checked this sensitivity with a test fit by increas- 
ing the errors on the measured dN/dF at low flux and 
our model did not prefer a different amplitude or shape 
to the source counts at the low flux where the efficiency 
for blazar detection is low. 

Our model fits the current DGRB, and, furthermore, 
predicts the DGRB for the expected enhanced sensitiv- 
ity to point sources after 5 years of Fermi-LAT data, 
2 x 1CT 9 photons cm~ 2 s _1 , which will resolve 94.7^2;?% 
of the flux from blazars. This expected enhanced point 
source sensitivity value is the Fermi-LAT Collaboration's 
estimate of the LAT flux sensitivity to point sources at 
high-latitude with gamma-ray index of ^2 29J . [30 In 
ABH2, we find the 68% and 95% confidence level (C.L.) 
upper and lower limit forecasts for the DGRB E 2 d<fr/dE 
when varying the fit parameters. This model finds that 
the DGRB will reduce by a factor of 1.6-2.6 (95% C.L.) 
with the spatial point-source resolution of the blazar con- 
tribution, after five years of the Fermi-LAT mission. The 
current and forecast spectra are shown in Fig. [T] In our 
forecasts, we calculate the limits using the flux as pre- 
dicted by the model at the minimal and maximal values, 
not simply performing a scaling of the limits. 

From the first to second Fermi source catalogs, there 
were 162 potentially spurious sources designated, indi- 
cating the sources' further identification with a spatially 
extended source, source variability, or other systematic 
effects [31 . The catalogs include sources with the test 
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FIG. 1. Shown are the current Fermi-LAT observed DGRB 
[3] in grey, and the forecast DGRB upper and lower 95% 
C.L. central values as the boxed (magenta) regions. Also 
shown are the expected emission from the Galactic and ex- 
tragalactic contribution for WIMP annihilation into bb for the 
WIMP particle masses 10 GeV and 100 GeV, for the canonical 
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statistic TS = 25, while the Fermi-LAT DGRB analysis 
only removed sources with TS = 50, or with higher signif- 
icance. (For the definition of TS, see Eq. 20 of Ref. [32].) 
This type of spurious contamination may alter forecasts 
for the DGRB, though the higher significance required 
for the exclusion of sources in the DGRB spectrum would 
likely reduce or eliminate this systematic effect. 



III. CURRENT AND FUTURE SENSITIVITY 
TO DARK MATTER MODELS OF THE 
FERMI-LAT DGRB 

The signal constrained by observations of the DGRB 
is the annihilation or decay of dark matter both in the 
Milky Way Galaxy and extragalactically. There is an 
irreducible contribution to the background from Galac- 
tic annihilation or decay that is isotropic and equal, at 
minimum, to the emission from the Anti-Galactic-Center 
(AGC). Here we examine in detail constraints on annihi- 
lating and decaying dark matter from both prompt and 
IC emission of photons. 



Diffuse Emission from Annihilating Dark 
Matter 



The products of dark matter annihilation emit in 
gamma-rays in several ways: annihilation channels that 
include the direct emission of a photon, decay of annihi- 
lation products into photons, and IC scattering of daugh- 
ter electrons off of background radiation. Through loop 
contributions to the annihilation cross section of dark 
matter, it is possible to have a direct 77 line signal, 
through a typically small branching fraction. If the dark 
matter annihilation particles include hadrons, then the 



decay chain of the products will lead to neutral pion de- 
cay into photons. Additionally, electrons among the dark 
matter annihilation daughters will up-scatter background 
photons from the cosmic microwave background (CMB) 
and starlight to gamma-ray energies through IC scatter- 
ing. We will consider the latter two cases in this paper: 
photons from prompt radiation in bremsstrahlung and 
hadronization, and IC emission from electron daughter 
particles. We calculate the photon and electron spec- 
trum from annihilation using the software P YTHI A [33] . 
To be conservative, we only consider IC emission from 
the CMB, not starlight or the infrared (IR) background, 
which is a good approximation in the direction of the 
AGC that we will consider. 

There are two sources of dark matter that contribute to 
the DGRB, Galactic and extragalactic dark matter. The 
signal from Galactic dark matter is largest in the line-of- 
sight toward the Galactic Center and is much smaller 
in other directions. However, there is an irreducible, 
isotropic signal from the Galactic dark matter that is 
equal to the signal from the AGC. This isotropic sig- 
nal can be the strongest dark matter contributor to the 
DGRB given conservative assumptions about the extra- 
galactic contribution [5]. 

To calculate the dark matter annihilation flux for a 
given cross section (<jav) and photon spectrum dN~/dE, 
we follow the treatment of Ref. [5] . The Galactic contri- 
bution to the diffuse flux is given by: 
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r ga i(b,£,x) = Jr\ - 2xR G cos{£) cos(6) + x 2 , (3.4) 

evaluated conservatively at the AGC (b = 0°, I = 
180°). In these equations, m x is the dark matter par- 
ticle mass, x is the line-of-sight distance, Rq is the dis- 
tance from the Galactic Center to the sun, and Jo = 
l/[8.5 kpc(0.3 GeV cm -3 ) 2 ] is an arbitrary constant that 
cancels in the final expression for flux. For the Fermi- 
LAT, the solid angle above \b\ > 10° has Att ohs = 10.4. 
The quantity J is the dark matter density squared inte- 
grated along the line of sight, and J'(AGC) is this value 
averaged over the observed angular sky region. For dark 
matter density p we use the minimal Einasto profile for 
the Milky Way halo: 



PEinastoM = Ps exp 
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(3.5) 



with ole = 0.22, r s = 21 kpc, r Q ~ 8.28 kpc, and 
Pq = 0.385 GeV cm~ 3 as in Ref. [6]. This profile is a 
conservative choice and gives J'(AGC) = 0.62, and ex- 
treme assumptions about the Milky Way dark matter 
profile only change this value by ~ 10%. 
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The Milky Way dark matter halo has abundant sub- 
structure which enhances the annihilation rate of dark 
matter. Following Ref. [31], the boost factor for the an- 
nihilation due to substructure is 
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The fraction f s of the halo volume is filled with a smooth 
dark matter component with density p%. The maximal 
density of the PDF p max is taken to be the scale den- 
sity p s of the earliest forming halos. The first term in 
Eq. (3.6), f s e A , is due to the variation in the smooth 



component, which contributes only a few percent to the 
boost, and therefore we ignore it. The second term is the 
boost factor due to substructure. The total luminosity 
boost due to the entire Galactic halo within radius R is 



B(< R) = 



§q B{r)p(r) 2 r 2 dr 
Jcf p( r ) 2r2 dr 



(3.7) 



where r is the halo-centric radial coordinate. The an- 
nihilation rate is larger in all directions than the AGC, 
therefore we calculate the luminosity with boost along 
that line of sight as the minimal annihilation rate due to 
our Galactic halo. Along the line-of-sight, 

.7boo S t(6, £) = Jo I S(r ga i(6, £, x))p 2 (r gal (6, t, x)) dx. 

(3.8) 

There is a partial reduction of the total luminosity boost 
[Eq. pTr) ] along the line of sight to the AGC due to 
our presence within the Galactic halo. Using the central 
value of a = from simulations, the boost is -Bagc = 
Jboost(AGC)/ l 7(AGC) = 3.3. Though the boost factor 
toward the Galactic Center is expected to be unity [34], 
that from the total Galactic halo can approach ~20 to 
2000. Therefore, the approximation of a total Galactic 
boost in the DGRB field of view as 3.3 is conservative. 

In addition to this Galactic contribution to the dark 
matter flux, we include a subdominant contribution from 
extragalactic dark matter annihilations [351 136] . This 
contribution is given by 
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where Hq is the Hubble constant, f2 m is the matter den- 
sity in units of the critical density, p cr ; t , and the fraction 
of matter in dark matter is /dm = ^dm/(^dm + &b) 1=3 
0.833, where we take the fraction of critical density of 
the dark matter as Odm = 0.237, and baryon density 
il h = 0.0456 [37]. Here, E' = E(l + z) is the source 
energy of the photons, and z up = (m x /E) — 1 is the 
maximum redshift to get a photon with energy E. The 



factor f(z) accounts for the increase in density squared 
during halo growth and the redshift evolution of the halo 
mass function. We adopt the fit of Refs. j3B][33], namely: 



f(z) = /o io°-9[-p(-o.9 Z )-i]- 



O.I62 



(3.10) 



For the Einasto profile, /o — 3 x 10 4 . We also include the 
boost factor of 6.6, from the total luminosity of a halo 
B(< R). This extragalactic contribution only accounts 
for < 30% of the total diffuse flux from dark matter. 

To calculate the limits on (crv), we attribute all the 
DGRB to dark matter annihilation (or decay), includ- 
ing both the Galactic and extragalactic dark matter con- 
tributions. This provides an upper limit on the cross 
section of annihilating dark matter (or a lower limit on 
the lifetime of decaying dark matter). We use the upper 
and lower 68% and 95% C.L. forecast fluxes correspond- 
ing to the extremal upper and lower fluxes in the three- 
dimensional contour in q, k and 71 parameter space, as 
constrained by the DGRB spectrum and source count 
distribution function, as described in our companion pa- 
per [14]. We take errors on the forecast DGRB for all of 
the upper and lower 68% and 95% C.L. forecast fluxes 
each to be proportional to its amplitude, which corre- 
sponds to the modeling methods of the DGRB measure- 
ment, though this is not necessarily the ultimate scaling 
of the errors. To do so, a Monte Carlo of the modeling 
methods of Ref. [3J would need to be performed, which 
is beyond the scope of this work. 

Figure [2] demonstrates the forecast for four canonical 
dark matter annihilation channels and how they compare 
to expected minimal supersymmetric model (MSSM) and 
minimal supergravity (mSUGRA) dark matter cross sec- 
tions. The cross sections are from a scan in MSSM and 
mSUGRA parameter space by Ref. [3jJ]. Note that our 
predicted constraints from Fermi-LAT's observations of 
the DGRB are comparable to current constraints Fermi- 
LAT observations of the Galactic Center [23], shown as 
a dashed line in Fig. [2] For comparison, in Fig. [2] we also 
show constraints from the Fermi-LAT analysis of stacked 
dwarf galaxies [40] [41] . The constraints shown are from 
Ref. [H]. In all relevant panels, we also show constraints 
from the HESS observation of the Galactic Center [12] as 
calculated in Ref. [4"3] . 

It has been proposed that a nonthermally-produced 
winolike dark matter annihilation could lead to the 
positron excess signal in PAMELA, which requires wino 
dark matter masses of 100 GeV < m x < 200 GeV [52]. 
As shown in Fig. [2](b) , these models are disfavored by 
the current constraints, and the forecast spectrum will 
be significantly more sensitive to this model. 

The standard thermal relic weakly-interacting dark 
matter annihilation cross section is ~ 3 x 10~ 26 cm 3 s _1 
[46] , which corresponds roughly with the broader 
MSSM/mSUGRA region. The plot of dark matter an- 
nihilating into t + t - , Fig. [2](d) , contains a region in pa- 
rameter space which has recently been claimed to be 
consistent with the morphology and spectrum of excess 
gamma-ray flux towards the Galactic Center [IS] , though 
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m x [GeV] m x [GeV] 

FIG. 2. Shown are our predictions for the Fermi-LAT sensitivity to constraints on dark matter in several canonical annihilation 
channels: (a) \X ~ ► bb; (b) XX ~ ► W + W~; (c) XX ~* tt; (d) XX t + t~ . The blue double and single hashed regions are 
the 68% and 95% C.L. predictions for 5-year Fermi-LAT sensitivity, respectively. Also shown is the limit from the current 
Fermi-LAT DGRB spectrum (solid blue). In panel (a), we also show the current constraints from Fermi-LAT observations of 
Draco (thick dotted line) [44 . In panels (a), (b) and (d), we show constraints from the stacking of dwarf galaxies (thin dotted 
line) [41 j - In the W + W~ channel, panel (b), the cross section versus mass for a nonthermal winolike neutralino is shown. The 
winolike dark matter of the PAMELA signal at m x ~ 200 GeV 22 is disfavored by the current constraints and will be further 
constrained in the forecast spectrum. The 'MSSM/mSUGRA' red-striped region is the expected cross sections for a sampling 
of points in supersymmetric parameter space |39| . The dashed line is a limit from the 3° x 3° in the Galactic Center |23| . In 
the t + t~ plot, panel (d), the green box is a region that could be consistent with an excess in the spectrum toward the Galactic 
Center 45 . In all panels, the triple-dot-dashed line is the limit from the HESS observation of the Galactic Center [35] for the 
respective channels in the case of a Navarro- Frenk- White halo profile, from Ref. [43] . 



such a signal is also consistent with emission from stel- 
lar clusters [37J . The width of the forecast region follows 
from the increasing width of the DGRB prediction with 
energy. Our model predicts a factor of 2 — 3 (95% C.L.) 
improvement in the sensitivity of the DGRB measure- 
ment to dark matter. 



B. Diffuse Emission from Annihilation of Dark 
Matter into Leptonic Modes and IC 

Dark matter that annihilates into leptons has been 
proposed as an explanation for an excess in cosmic-ray 
positrons seen by the PAMELA Collaboration [50] and 
a feature in the cosmic-ray e + /e~ spectrum from the 
Fermi-LAT [5T]- The DGRB is forecast to be sensitive 



to direct production models that fit the PAMELA 

excess and Fermi-LAT e + /e~ feature, such as Asymmet- 
ric Dark Matter [21]. Such direct annihilation models 
are already strongly disfavored by constraints of obser- 
vations by the Galactic Ridge by HESS [SI [SU |S3] , obser- 
vations by Fermi-LAT toward the Galactic Center |23| . 
and HESS toward the Galactic Center 021 [43]. The latter 
two constraints shown in Fig. [3] 

It has also been suggested that Sommerfeld-enhanced 
dark matter annihilation into four leptons, via a light 
mediator particle, could explain these signals. There- 
fore, we also consider how an improved DGRB will affect 
these leptonic dark matter annihilation channels. Note 
that such models are also constrained by detailed calcu- 
lations of the relic abundance from dark matter produc- 
tion in the early Universe and halo shapes [53], as well 
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FIG. 3. Shown in panel (a) are our predictions for Fermi-LAT sensitivity to n + channel dark matter annihilation. The dot- 
dashed line is the 95% C.L. limit on prompt and IC emission from Ursa Minor [44 . In this panel, we show the constraint from 
the stacking of dwarf galaxies (thin dotted line) [41]. In panel (b), dark matter annihilation into four muons via intermediate 0.3 
GeV scalar particles (j>. The blue double and single hashed regions are the 68% and 95% C.L. predictions for 5-year Fermi-LAT 
sensitivity, respectively. Also included is the limit from the current Fermi-LAT DGRB spectrum (striped blue region). The 
solid black line shows where the exclusion would be without the IC contribution. The light pink shaded region is consistent 
with a dark matter interpretation of the PAMELA signal and the dark red shaded region is consistent with a dark matter 
interpretation of the Fermi-LAT e + /e~ feature [481 149] . The triple-dot-dashed line is the limit from the HESS observation of 
the Galactic Center [35] for the respective channels in the case of a Navarro- Frenk- White halo profile, from Ref. |43| . In panel 
(b), the dotted line is a limit on dark matter annihilation from radio synchrotron from the Galactic Center (48]. The dashed 
line is the 99% C.L. limit from the 3° x 3° region toward the Galactic Center |23| . 



as distortions of the CMB spectrum [55]. Recent work 
finds that several reasonable parameter choices in mod- 
els of Sommerfeld-enhanced dark matter for PAMELA 
and the Fermi-LAT e + /e~ spectral feature avoid these 
limits [56| . These Sommerfeld-enhanced models are par- 
tially constrained by the current limits and will be further 
constrained by our forecast limits, as shown in Fig. J3jb). 

Dark matter annihilation into leptons produces fewer 
photons than the quark or gauge boson channels. How- 
ever, such annihilations do produce highly-boosted elec- 
trons which undergo IC scattering on the CMB and 
starlight which is then observable in gamma-rays. We 
calculate this IC contribution as in Ref. [23]. To be con- 
servative, we only include the scattering from the CMB, 
not from starlight or the IR background. 

The spectrum of IC photons coming from one dark 
matter annihilation is given by [23] [2Z] 



dN 
dE 



1 
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(3.11) 



where V(E, e) is the differential power emitted into pho- 
tons of energy E by an electron with energy e, £ denotes 
the total rate of electron energy loss due to IC scatter- 
ing, and Y(e) is the number of electrons generated with 
energy larger than e in one annihilation. To get the en- 
ergy of the annihilation products, we use the software 
PYTHIA 33J. In the Thomson limit, which is a very 



good approximation for CMB photons, one obtains 
4 
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(3.12) 



viE > e)= w E L dy 



x (2ylny + y + l-2 2 / 2 ), (3.13) 

where 7 is the Lorentz factor of the electron, e' is the 
energy of the initial CMB photon, y = E/(4'j 2 e / ), ot — 
0.665 barn is the Thomson cross section, and the radia- 
tion density of CMB photons at Tomb — 2.725 K is given 
by 



^CMB(e') 



1 



7T 2 exp(e'/TcMB) - 1 



(3.14) 



The flux coming from the IC contribution can then be 



calculated as in Eqs. (3.1 1-< 3.10 ) 



It is important to note that we are neglecting the diffu- 
sion of the electrons and positrons from the dark matter 
annihilation to the point where IC scattering occurs. As 
stated in Ref. |23| . this is a good approximation in re- 
gions away from the Galactic Center, as we consider in 
this work. For an analysis where diffusion is included, 
see Ref. [33] • To be more explicit, in Ref. [57], it was 
shown that the electron/positron flux produced by dark 
matter annihilations was only slightly modified by diffu- 
sion effects for distances from the Galactic Center larger 
than 8 kpc. Corrections were found to reach a factor of 
2 only at the lowest energies around 100 MeV (see their 
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Fig. 14). Therefore, diffusion effects are not significant 
in our calculation of the IC component. Note that the IC 
component is modeled in the measurement if the DGRB 
independent of the observation [3J, so there is no acci- 
dental subtraction of a potential IC signal. 

Our calculations for the standard p + pT lcptonic chan- 
nel are shown in Fig. [3]Ja). Also shown are the re- 
gions which are consistent with the PAMELA excess and 
Fermi-LAT e + /e~ spectral feature from Ref. [35], modi- 
fied for a higher local dark matter density of our minimal 
Einasto halo model (p@ — 0.385 GeV cm -3 ), and the 
local boost of B(r = R Q ) = 1.57 [Eq. ( pLoj )] from sub- 
structure. For this boost, we use the same parameters for 
the subhalo PDF as our subhalo boost annihilation sig- 
nal. Such models are already highly constrained by sev- 
eral gamma-ray observations as shown in Fig. J3ja) , and 
such models will be further constrained with our forecast 
DGRB sensitivity. Note that the IC contribution im- 
proves the bounds by several orders of magnitude. The 
IC gamma-ray flux contribution peaks at much lower en- 
ergies than the prompt component for a given dark mat- 
ter particle mass, which is why the width of the forecast 
region decreases when the IC component becomes impor- 
tant. Also shown in Fig. [3] are complementary limits on 
the cross section of /Lt + /x _ -channel annihilating dark mat- 
ter from other work. Shown in Fig. [3jb) is dark matter 
annihilation into four muons via intermediate scalars <f> 
with 0.3 GeV masses. Even in this less-constrained case, 
our forecast is that the Fermi-LAT measurement of the 
DGRB will have the sensitivity to detect or rule out a 
portion of the parameter space the dark matter interpre- 
tations of PAMELA. 



C. Diffuse Emission from Decaying Dark Matter 

Decaying dark matter is another, less constrained, pos- 
sibility for the source of the positron fraction signal seen 
in PAMELA and the Fermi-LAT e+ je~ spectral fea- 
ture [23] . To calculate the flux from decaying dark matter 
with lifetime r and photon spectrum dN^/dE, the proce- 
dure is very simil ar to the annihilating dark matter case. 



However, in Eq. (3.3|, p — > p and there is no boost fac- 



tor due to the lack of decay enhancement with density. 
Here, Eq. (3.1| becomes 



d<f>-y _ 1 Jagc 1 dN^ 
dEdQ t Jq 47tto x dE 



(3.15) 



Similarly, to calculate the extragalactic flux from decay- 
ing dark matter, in Eq. (3.9) replace 



(0"A«) (/DM^m) 2 Pcr 



1 (/DM^m)Pcrit 

t m Y 



(3.16) 



and drop all boost factors. 

Fig. [J] shows how the predicted DGRB value will con- 
strain the lifetime of a dark matter particle which de- 
cays into p + pT . The improved constraint will have the 




1000 



m x [GeV] 



FIG. 4. Our predictions for Fermi-LAT sensitivity to dark 
matter decaying into p + p~. The dark and light blue hashed 
regions are the 68% and 95% C.L. predictions for 5-year 
Fermi-LAT sensitivity, respectively. Also included is the limit 
from the current Fermi-LAT DGRB spectrum (blue striped 
region) . The solid black line shows where the exclusion would 
be without the IC contribution. The light pink shaded re- 
gion is consistent with a dark matter interpretation of the 
PAMELA signal, and the dark red shaded regions are con- 
sistent with a dark matter interpretation of the Fermi-LAT 
e + /e~ feature at 3- and 5-er [49] . The dashed line is from 
constraints on prompt and IC emission from dark matter an- 
nihilation made from Fermi-LAT observations of the Fornax 
cluster of galaxies . 



sensitivity to exclude or detect the decaying dark mat- 
ter interpretation of the Fermi-LAT e + je~ feature and 
should provide strong limits on an interpretation of the 
PAMELA excess gSISU]. The DGRB limit is compara- 
ble to and is forecast to be more sensitive than the limits 
on decaying dark matter from Fermi-LAT observations 
of clusters of galaxies [S5J , as shown in Fig. [4] 



D. Comparison to Direct Dark Matter Detection 
Limits on Light Dark Matter 

Dark matter may be detected through two distinct 
methods: indirect-detection experiments seek the anni- 
hilation or decay products from dark matter in the Uni- 
verse, while direct-detection experiments look for the re- 
coil of heavy nuclei after their collision with a dark matter 
particle from our Galactic halo. In general, the interac- 
tion cross section between annihilating dark matter is 
not simply related to the interaction cross section be- 
tween dark matter and nucleons. A few recent direct- 
detection experiments, however, have seen signals that 
could be caused by a light dark matter particle, includ- 
ing DAMA [BIMIEI], CoGeNT [60] and CDMS [64]. 

These light dark matter signals could be due to dark 
matter that interacts through the exchange of Higgs 
bosons [53K2S]. For such a dark matter candidate, the in- 
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m x [GeV] 



FIG. 5. Our predictions for Fermi-LAT sensitivity to light 
dark matter coupled via Higgs-like couplings |25| . The dark 
and light purple hashed regions are the 68% and 95% C.L. pre- 
dictions for 5-year Fermi-LAT sensitivity, respectively. Also 
included is the limit from the current Fermi-LAT DGRB spec- 
trum (black line). To the right of the dotted pink line is 
exclusion from the XENON100 Collaboration [52]. The red 
"CoGeNT" region is consistent with the findings of the Co- 
GeNT Collaboration [60]. The blue upper "DAM A" region 
is the DAMA signal without channeling [61]. The orange 
lower "DAMA*" region is the DAMA signal if channeling is 
included 62, 63 . The region labeled "CDMS" is excluded by 
that experiment's light dark matter search at 95% C.L. |64| . 



direct annihilation cross section and direct nuclear cross 
section are related by 

- v A? m 2 f(m% - ml) 3 / 2 

<r ind (SS^ff)-=n c ^ /V S 4 3 /; (3.17) 



<T dir (SiV >SN) A ' f. 



2™2 



N 



(3.18) 



b h" b S 



^ Cdir c ^ J 2 m z N Hpm s 



where m/j is the Higgs mass, A^ is the dark matter-Higgs 
coupling, n c = 3(1) for quarks (leptons), fj, r is the dark 
matter-nucleon reduced mass, / ~ 0.3, and ms is the 
dark matter mass [25 . The sum in Eq. (3.19) is over 



all annihilation products, which are dominated by the 
6-quark, the c-quark, and the r-lepton. Through this 
ratio, we can relate our predicted DGRB limits on indi- 
rect detection into limits on direct-detection experiments. 
For such a Higgs-mediated dark matter model, we com- 
pare our projected limit to the findings of several direct- 
detection experiments in Fig. [5] 

The limit on the dark matter-nucleon cross section as 
found by our DGRB forecast is competitive with the lim- 
its by the XENON100 Collaboration in the lowest mass 
range [59 . The current Fermi-LAT DGRB values rule 



out the DAMA region without channeling and some of 
the region consistent with a dark matter interpretation 
of CoGeNT. After a Fermi-LAT 5-year run, we forecast 
the DGRB spectrum to have the sensitivity to exclude 
most of the CoGeNT region consistent with dark mat- 
ter interpretations in the spin-independent case. This is 
complementary to the findings of direct-detection exper- 
iments since the DGRB indirect-detection limits tend to 
exclude lower dark matter masses than direct-detection 
experiments. 



IV. CONCLUSIONS 

The DGRB as measured by Fermi-LAT is one of the 
most powerful constraints on annihilating weak-scale par- 
ticle dark matter. We show that the likely resolution 
of blazars into point sources by Fermi-LAT — and their 
automatic removal from the DGRB measurement — will 
enhance the sensitivity of the DGRB to dark matter an- 
nihilation by a factor of 2 to 3 (95% C.L.), depending on 
the channel, mass scale, and true realization of the blazar 
distribution and SED sequence. 

We find the forecast dark matter sensitivity of the 
DGRB observation to both prompt and inverse-Compton 
photon emission to be comparable in sensitivity with 
other limits on annihilating weak-scale dark matter. The 
DGRB is forecast to be comparable to current limits from 
Fermi-LAT observations toward the Galactic Center [33] , 
individual dwarf galaxies [44], and, in the case of decay- 
ing dark matter, observations of clusters of galaxies [55] . 
This sensitivity makes the DGRB among the best meth- 
ods of detecting or constraining dark matter with the 
Fermi-LAT mission. The forecast for the DGRB projects 
it to be less constraining, for certain channels and particle 
masses, than current stacked dwarf analyses with Fermi- 
LAT [40, 41j and observations of the Galactic Center by 
HESS [HiS]. 

The future resolution and reduction of the DGRB into 
blazar point sources highlights and enhances the possi- 
bility that the Fermi-LAT experiment will either detect 
or constrain the dark matter in a robust yet conservative 
way. 
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